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ABSTRACT 



The nature of the binary systems giving rise to Type la supernovae (SNela) 
remains an unsolved problem. In this Letter we calculate, from the statistics 
of initial conditions (masses and binary separations), the mass, luminosity, 
and velocity distributions of the possible binary companions (main-sequence 
star, subgiant, red giant) following the explosion of the white dwarf which 
gives rise to the SNela. Those companions could be detected from either their 
proper or their radial motions, by means of high-precision astrometric and 
radial-velocity measurements in young, nearby supernova remnants. Peculiar 
velocities typically ranging from 100 to 450 km s _1 should be expected, which 
places proper-motion measurements within reach of HST instruments and makes 
radial-velocity ones feasible with 2.5-4m class telescopes from the ground. 
Detections would solve the long-standing problem of which kind of binaries 
do produce SNela and clear up the way to accurate physical modeling of the 
explosions. 



Subject headings: supernovae: general — binaries: close — ISM: supernova 
remnants — astrometry 



- 3- 



1. Introduction 

Type la supernovae (SNela) are unanimously attributed, nowadays, to the 
thermonuclear explosion of a white dwarf (WD) due to accretion of matter from a 
companion in a close binary system. The mass-accreting WD is also generally thought to 
be made of carbon and oxygen (CO WD). There the agreement stops, and different kinds 
of systems are proposed as candidates to SNela progenitors (see Iben 1997, for a review). 

The nature of the WD companion, the mass-donor, is unclear: it could either be 
another WD (double-degenerate or DD systems) or a still thermonuclearly evolving star 
(single-degenerate or SD systems). The evolutionary stage of the companion in the SD 
systems could be anything from a main-sequence star to a supergiant, and mass accretion 
proceed via Roche-lobe overflow (driven by thermonuclear evolution of the companion, 
magnetic braking, or gravitational wave radiation) or by capture of matter from stellar 
wind. All that also bears on the explosion mechanism: on how the explosion is triggered and 
how it develops. Those uncertainties still raise doubts as to the use of SNela as calibrated 
candles to probe the dynamics of the Universe, the calibration being purely empirical and 
based on the nearby SNela sample (see Drell, Loredo, & Wasserman 2000, for instance). 

Recently, evidence has been gathered against the merging of WDs giving rise to SNela, 
from the high efficiencies found in producing SNela out of star formation. These efficiencies 
have now been measured at several redshifts and appear to be similar, every ~ 700 M e 
going into star formation giving rise to 1 SNela (Ruiz-Lapuente & Canal 2000). 

No previously observed system has ever been identified as the progenitor of any 
SNela. The historical SNela took place at a time when astronomical instrumentation was 
rudimentary. They can still give us, nonetheless, important and maybe decisive clues on the 
nature of the long-sought SNela progenitor systems (Ruiz-Lapuente 1997). If they were 
DD systems, there should be no companion left since it is destroyed in the accretion process 
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already. On the contrary, in SD systems the companion can survive the explosion in most 
cases and show characteristics allowing its identification. 

Of the historical SN, from the records on lights curves only two have been identified as 
SNela: SN 1006 and Tycho's SN (SN 1572). From their X-ray morphology (spherically- 
symmetric shell), X-ray spectra (showing high abundances of Fe, Ni, Si, Ca, Ar), Fe 
enrichment in their optical filaments, and lack of any central X-ray source, a number of 
nearby, young Galactic supernova remnants (SNR) can be attributed to SNela (Ruiz- 
Lapuente et al. 2000). As we will see, precise enough astrometric and radial-velocity 
measurements of the central regions of such SNR should allow detection of the SNela 
companions from either their large proper motions or high radial velocities, due to 
disruption of the binary orbit plus kick from the impact of the SN ejecta on them. For SNR 
ages <, 10 4 yr, the companions should still exhibit in their outer layers the consequences of 
the strong perturbation of thermal equilibrium from partial stripping of the envelopes and 
deposition of energy from the ejecta into the layers that remain bound, and they may also 
show chemical contamination from the SNela ejecta. 

2. Modeling, Results, and Discussion 

The SD scenarios thus far proposed to explain the origin of SNela involve different 
kinds of systems. We will only consider Ch-mass explosions, which are generally thought 
to produce the bulk of SNela, so the SN ejecta impinging on the companion will always 
be identical (differences due to the explosion mechanism should be minor and will be left 
aside). Depending on the kind of system, however, the mass and evolutionary stage of the 
companion, as well as the separation between the two stars at the time of explosion, will 
be different and that will determine the space velocity and the mass, luminosity, surface 
temperature, and surface abundances of the surviving star. 
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We will consider the cases of a main-sequence companion filling its Roche lobe due 
to angular-momentum loss via magnetic braking (cataclysmic variable or CV system), 
and that of a subgiant or red-giant companion filling its Roche lobe due to thermonuclear 
evolution (cataclysmic-like or CLS system). The latter are thought to show up as luminous 
supersoft X-ray sources before exploding as SNela. In the modeling of the CLS systems, 
the stabilizing effects on mass transfer of the strong wind generated in the WD surface have 
been taken into account as in Hachisu, Kato, and Nomoto (1996). The cosmic SNela rates 
thus predicted give the best agreement with the observations (Ruiz-Lapuente & Canal 
2000). 

Our purpose here is to derive the statistical distribution of velocities, masses, and 
luminosities of the companions after explosion. We first run and updated version of the 
Monte Carlo scenario code we used to model SNela rates in previous work (Ruiz-Lapuente, 
Burkert, & Canal 1995; Canal, Ruiz-Lapuente, & Burkert 1996; Ruiz-Lapuente & Canal 
1998), to obtain the distributions of companion masses and orbital separations at the time 
of explosion, for the systems considered. Those are displayed in Figure 1 for the CV systems 
(main-sequence, MS companion) and the CLS systems (distinguishing between subgiant 
SG, and red-giant RG companions). Note, in the upper panel, that the mass distributions 
of the SG and RG stars in the CLS scenario are the same, the SG case corresponding to 
the narrow range of final separations shown in the lower panel whilst the RG case covers a 
much broader range, with an almost flat distribution (same panel). In the CV scenario, 
the much narrower distribution of final masses, peaking at ~ 0.6M Q , means an also narrow 
distribution of final radii, the companion being still in the MS. Making contact with the 
Roche lobe while having a companion (the WD) that has always the Ch-mass produces the 
sharp peak in the distribution of final separations that we see in the lower panel. 

From the initial conditions above, the effects of the impact on the companion are 
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calculated in a semianalytical approximation equivalent to that used by Wheeler, Lecar, & 
McKee (1975), based on previous work by Colgate (1970). The companion star is divided 
into concentric cylinders and the ejecta are treated as a plane slab. If the momemtum 
incident on a cylinder can accelerate it to the scape velocity, the corresponding mass 
is directly stripped. A strong shock wave is driven into the rest of the star: that will 
impart momentum to it (kick) and the heating will make the internal energy to exceed the 
gravitational potential energy in extra layers which will thus be expelled. The approximation 
is made that the fluid velocity just behind the shock is equal to the mean velocity of all 
the shocked stellar material. The energy per unit mass rapidly increases as the mass per 
unit area in the cylindrical shells decreases: energy deposition is thus concentrated in the 
outer layers, which ensures that the companion is never entirely disrupted. The results have 
been compared with the two-dimensional hydrodynamic simulations of Marietta, Burrows, 
& Fryxell (2000). Especifically, this has been done for the cases of a 1 M Q MS companion 
at separation a = 3 f?©, a 1.1 M & SG at a = 4.9 R &) a 2.1 M Q SG at a = 6.5 R e , and 
a 1 M RG at a = 430 R®. The agrement (maximum discrepancy below 20%) is good 
enough to confirm our approximation as a reliable tool for the statistical evaluation of the 
characteristics of the SNela companions after explosion, which would be hard to do by 
means of detailed hydrodynamic calculations. Note that the combinations of companion 
masses and separations at the time of explosion in Marietta, Burrows, & Fryxell (2000) 
(see their Table 4) are just adapted from examples proposed by Livio & Truran (1992) (CV 
system) and Li & van den Heuvel (1997) (CLS system), whereas in the present work the 
actual distributions of those parameters have been calculated for each scenario. That is why 
the cases above (with the exception of the 1.1 M & SG companion at a = 4.9 R Q ) are not 
typical of the predicted distributions. However, Marietta, Burrows, & Fryxell (2000) have 
compared they results with an aproximation similar to ours and find best agreement for 
binaries near Roche lobe overflow, which is always our case. That reinforces the significance 
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of our comparisons. 

The structures of the MS and SG companions prior to explosion have been obtained 
by homologous transformations of a 1M Q model in the corresponding stages. Typically, 
10-20 % of the companion mass is lost. The RG companions are modelled as a compact 
He core surrounded by an extended H-rich envelope, the latter having the structure of a 
n — 1.5 polytrope. The relationship between RG mass, radius at the time of explosion, 
and core mass, has been taken from the stellar evolution calculations by Politano (1988), 
which cover the initial mass range 0.1 M < M < 10 M & . We find that in all cases the RG 
envelope is completely stripped off, as in the 2D simulations of Livne, Tuchman, & Wheeler 
(1992). Marietta, Burrows, & Fryxell (2000), in their recent 2D hydrodynamic modeling 
find, however, that a small fraction of the envelope is retained, which would greatly help 
in the identification of possible SNela companions when exploring the central regions of 
young, nearby SNR, as discussed below. 

The kick velocities imparted to the companions by collision wih the SNela debris are 
also calculated, but they are always much smaller than the typical orbital velocities at 
the time of explosion, for the three types of companions (MS, SG, RG). The predicted 
distributions of total velocities (orbital plus kick) are shown in Figure 2. 

The results corresponding to the most typical combinations of mass and orbital 
separation, for each type of companion, are shown in Table 1. In column 2 we give the 
apparent visual magnitude m v expected (for distances of 1 kpc and 5 kpc). In column 4, 
the average radial velocities v r (equal to the average transversal velocities v t , since the total 
velocities should be distributed at random with respect to the observer's line of sight). In 
column 2 is the proper motion fi (in arcsec yr _1 ), and in column 5 the maximum angular 
distance 9 (in arcmin) from the explosion site that would be reached 10 3 yr after the SN 
event. Given in the footnotes to Table 1 are the pre-explosion and post-explosion masses 
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for the three cases. We must stress here that the apparent visual magnitudes given in Table 
1 are calculated from the luminosities that the companions would reach after relaxing to 
thermal equilibrium, and thus are lower limits to the actual brightnesses to be expected up 
to a few thousand years after explosion. It is clear that the strong shock wave will not only 
unbind a fraction of the outer layers and give a global kick to the companion but it will 
also strongly perturb the hydrostatic and thermal equilibrium of the fraction that remains 
bound. Mechanical equilibrium should be restored on a hydrodynamic time scale and 
thermal equilibrium on a longer, Kelvin-Helmholtz time scale. Indeed, Marietta, Burrows, 
& Fryxell (2000) find in their 2D numerical simulations that after impact the MS and SG 
companions would swell up, their luminosities rapidly increasing by factors of up to a few 
thousand, to relax back to luminosities of the order shown in Table 1 in ~ 10 4 yr. The RG 
companions would retain a residual H-rich envelope, still able to feed the H-burning shell 
and would evolve away from the RG branch at constant luminosity and increasing effective 
temperature, on a time scale of ~ 10 4 — 10 5 yr. Therefore, SNela companions left by recent, 
nearby explosions, should in fact be detectable at magnitudes significantly brighter than 
those shown in Table 1. To illustrate this we also give, within parentheses, the magnitudes 
expected for the MS companions (the less luminous ones) at t ~ 10 3 yr after explosion. The 
peculiar velocities predicted are one order of magnitude larger that the systemic velocities 
of field stars at the corresponding distances. The angular sizes of the regions to be explored 
(given by 9 in Table 1) ensure that the number of stars bright enough to be candidates to 
SNela companions, contained within their boundaries, is quite small. Detection of proper 
motion at the level of 0.01 - 0.005 arcsec yr -1 , even at the magnitudes estimated in Table 
1, is within the reach of the instruments FGSlR(brighter targets) and WFPC2 (fainter 
targets) on board the Hubble Space Telescope (Ruiz-Lapuente et al. 2000). Failure to detect 
any candidates would be strong evidence that none is left after the SNela explosions and 
almost completely discard the SD system hypothesis. 
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The possible candidates selected from their high proper (or radial) motion should then 
be spectroscopically analyzed. As we have seen, MS and SG candidates should be swollen 
and overluminous and RG candidates would have high effective temperatures. Moreover, 
Marietta, Burrows, & Fryxell (2000) find that whereas no material from the ejecta should 
be accreted during impact, there is a possibility of later accretion of material enriched in Fe 
and intermediate-mass elements from the low-velocity tail of the ejecta, long after impact. 

There are still many uncertainties concerning close binary evolution, binary mass 
transfer, and the process of mass growth of a WD from accretion until it reaches the 
Ch-mass. To which extent could those uncertainties affect our result that any SNela 
companion should be moving at high velocities after the explosion and be luminous enough 
to be detected in the central regions of young, nearby SNR? 

In the case of the CV systems, the companions are MS stars that continuosly fill 
their Roche lobes due to angular momentum loss by magnetic braking (plus emission of 
gravitational waves when the orbits become very close). The WD mass being equal to the 
Ch-mass, separation is only a function of companion mass, which in turn is just the initial 
mass of the star minus the fraction transferred to the WD and the (minor) fraction lost to 
stellar wind. The efficiency of magnetic braking is here the main unknown, but it hardly 
affects the distribution of final masses (and thus separations) shown in Figure 1 nor the 
values predicted in Table 1. 

Although the CLS systems involve companions in a more advanced evolutionary stage 
than the CV systems, our main conclusions are equally robust. When the companion is 
a SG, the case is as clear-cut as for the CV systems. For the RG companions, the only 
worrying possibility, that the companion velocities actually were much lower than in our 
modeling, can be reliably excluded: since v or b oc dJ 1 ^ 2 (Mi + M 2 ) 1 ^ 2 and the WD mass 
is always equal to the Ch-mass, we approximately have v orb oc a7 1 ^ 2 (M 1 + 1.4) 1 / 2 , Mi 
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being the companion's mass. Increasing ay (and, less efficiently, lowering Mi) would thus 
decrease v or b- But the companion must be filling its Roche lobe Rl, and we also have 
(Paczyhski 1971) that Rl oc a/ (0.351 + log Mi) (where in the numerical constant we have 
again taken M 2 ~ M C h)- Increasing a/ should, therefore, increase Rl faster than it would 
decrease v or b (decreasing M 1 cannot go very far, since the companion's initial mass must 
be Mi >, 0.8 M & in order to have evolved from the MS in less than a Hubble time), and 
substantial increase postpones Roche-lobe overflow to the red supergiant (RSG) stage, 
when a common envelope would form and be ejected, that suppressing mass transfer to the 
WD altogether and leaving a DD system, no longer a SD one. 

A class of SD system that we have not considered here is the one in which the WD 
accretes H-rich material from the wind of a RG or RSG that is not filling its Roche lobe, 
a symbiotic system (SS). There the orbital separation could be larger and the orbital 
velocities smaller. However, the difficulty in making the WD grow to the Ch-mass through 
such a low-efficiency accretion process, avoiding explosive ignition of either H or He in the 
outer layers (which requires high accretion rates), has discarded SS systems as possible 
progenitors of the bulk of the SNela at least (Kenyon et al. 1993). 

Finally, that the RG companions do lose their envelopes at most, and that the SG 
and MS companions typically just lose a minor fraction of their masses upon being hit by 
the SNela ejecta, is consistently found in all hydrodynamic simulations and semianalytical 
approximations. That upon impact and stripping of the outer layers the new surface layers 
of the companions will be out of thermal equilibrium, and that the objects should thus 
initially appear overluminous as compared with the equilibrium state they will approach on 
a thermal time scale, seems hardly disputable as well. 



3. Conclusions 



We have added to the results from existing hydrodynamical studies of the impact of 
SNela into their companion stars the statistical evaluation of the distribution of velocities to 
be expected for different kinds of companions. We have shown that they should be moving 
at velocities much higher than those of the surrounding stars: for RG companions there is a 
probability P > 90% that the velocities are higher than 100 km s' 1 ; for SG P > 99% that 
v > 200 km s^ 1 , and for MS P > 99% that v > 450 km Such velocities make their 
identification possible in the central regions of young, nearby SNR of SNela origin. For 
distances d < 5 — 8 kpc and SNR ages <, 10 3 — 10 4 yr, detection of proper motion is within 
the reach of instruments on board the HST, and that of radial motions can be made with 
2.5-4m class telescopes from the ground. The angular sizes of the regions to be surveyed 
are small enough that the numbers of stars within them, bright enough to be candidates to 
SNela companions, remains low. There are 8 Galactic SNR (including those of SN 1006 
and SN 1572) that most likely formed from SNela explosions and whose distances and ages 
allow thorough searches (Ruiz-Lapuente et al. 2000). 

Candidates selected from their peculiar velocities and positions close to the centers of 
the SNR could be confirmed by spectroscopic characteristics arising from their being out 
of thermal equilibrium due to heating by the impact of the SN debris, according to recent 
hydrodynamic simulations (Marietta, Burrows, & Fryxell 2000). Contamination by material 
accreted from the low-velocity tail of the ejecta is also possible. 

Our conclusions are largely independent from the unknowns still affecting close binary 
evolution, mass transfer, and growth of WDs up to the Ch-mass. 

Therefore, it now appears feasible to test the SD system hypothesis for the origin of 
SNela by means of observations of Galactic SNR. Detection of confirmed companions would 
not only validate the SD system hypothesis, but it would also allow physical modeling of 
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SNela on a firmer basis and, through that, make much more reliable than at present the 
use of SNela in Cosmology. 
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Table 1: Typical apparent magnitudes, proper motions, radial velocities, and maximum 
angular distance from explosion site (after 10 3 yr) of SNela companions 

Companion type my /i (arcsec yr -1 ) v r (km s _1 ) 6 (arcmin) 

and distance 



Main sequence a 



1 kpc 18.7 (11.2)* 0.067 320 1.6 

5 kpc 22.2 (14.7)* 0.013 320 0.3 

Subgiant b 

1 kpc 12.6 0.038 180 0.9 

5 kpc 16.1 0.008 180 0.2 

Red giant c 

1 kpc 10.5 0.015 70 0.4 

5 kpc 14.0 0.003 70 0.1 



a Pre-explosion mass: 0.60 M ; post-explosion mass: 0.51 M . b Pre-explosion mass: 
1.00 M ; post-explosion mass: 0.85 M . c Pre-explosion mass: 1.00 M ; post-explosion 
mass: 0.33 M Q (electron-degenerate He core). 
* Expected at t ~ 10 3 yr after explosion. 
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Fig. 1. — Predicted distributions of SNela companion masses (upper panel) and of orbital 
separations (lower panel), just prior to SNela explosion, for the two types of single-degenerate 
systems considered: a WD in a cataclysmic variable system, where the companion is a MS 
star (dotted line), and in a cataclysmic-like system, where the companion can be either a SG 
(dashed line) or a RG star (solid line). Note that in the lower panel the values corresponding 
to the RG case have been multiplied by a factor 5, to improve legibility. 
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Fig. 2. — Predicted velocity distributions of possible binary companions of the exploding 
WD, following the SNela explosion, for the same single-degenerate systems as in Figure 1. 
Our statistical evaluation suggests that if the companion star is a RG there is a probability 
P > 90% that it will move at a velocity larger than 100 km s" 1 . For the SG, there is a 
P > 99% that it will move at a velocity larger than 200 km s -1 , and for the MS companion 
there is probability P > 99% that it will be moving at a velocity larger than 450 km s _1 . 
The chances that the central regions of young SNR produced by SNela contain a star moving 
at a velocity larger than 100 km s _1 are thus very high if the SD scenario is true. 



